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It is shown that the BKZ incompressible elastic fluid theory is in excellent agreement with experi- 
mental results obtained in simple extension. From single step stress-relaxation data, the stress-strain 
response for a number of other simple extension histories are calculated from the theory and are 
compared with experiments. 
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1. Introduction 

In 1963 Bernstein, Kearsley, and Zapas [l] 1 pre- 
sented a theory of an elastic fluid which gave the 
correct stress-relaxation response for a large variety 
of elastomeric materials, including vulcanized rub- 
bers. A principle attractiveness of this theory is 
its relative simplicity; with a single integral in time, 
it describes the stress-strain behavior for all types 
of deformation histories. In the case of simple ex- 
tension; it predicts the behavior in any uniaxial strain 
history from the results of single step stress-relaxation 
experiments which cover the same range of extension 
and time. We designed a series of experiments to 
check the validity of this theory and found, as is 
shown in section 4 of this paper, excellent agreement 
with experiment in all cases. We are aware that ex- 
periments cannot prove a theory. From our results, 
however, we feel strongly that a single integral expres- 
sion with a nonlinear integrand such as the BKZ elastic 
fluid equation is sufficient to describe the stress-strain 
behavior of elastomeric materials. 



2. General Considerations 

For an extensive description of the BKZ elastic 
fluid, we refer the reader to the initial papers [1, 2J. 

Here we shall use only the definitions and deriva- 
tions necessary for the description of uniaxial defor- 
mations, i.e., simple extension of an incompressible 
material. 

A particular motion of the material may be specified 
by a relation between the Cartesian components of 
the position of a particle at time t, Xt(t) and of its posi- 
tion at time r, -*";(t). 



Xi(t) = Xi(xk(r), t, t) 



1 Figures in brackets indicate tlu 



i. A = 1. 2. 3. 
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The relative deformation gradients Xik(t 9 r) are defined 
by 



Xik(t, t)-- 



dxjt) 
Bxkir) 



The left Cauchy-Green tensor dj{t, r) is defined by 

dj(t, T) = Xik(t, T)Xjk(t, t) 

where, in accordance with Cartesian tensor notation, 
repeated indices indicate summation over the values 
1, 2, 3. The principal invariants of q/U, r) are 

h{t, T) = |*+fi+£S 

lit, T) =&!+&* + &* 

h{ t ,T)=e,m 

Ai(«) 



, , Aito , 
where £1 = —7 -, £> 
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In view of the condition of incompressibility which 
requires that I 3 (t, r)= 1, the constitutive equation 
for the BKZ elastic fluid can be written 



n j = -p8, + 2j_^-c ij tor)--c^tor) 



dr (1) 



where dj are the components of the stress tensor, p is 
some hydrostatic pressure, U is a function of /i, I>, 
and t — r, 

U=U(I 1 (t,T)J 2 (t,T\t-T) 

and cj. l {t, t) are the components of the inverse of c^tj). 
We may describe isochoric uniaxial strain in a co- 
ordinate system in which the %\ direction corresponds 
to the direction of stretch. Then the particle which 
in some reference configuration is located at Xi, X2, X$ 
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has a position at time t, x(t) given by 

Xi(t)=k(t)X 1 

x 2 (t) = k-^(t)X 2 

x 3 (t)= X-^Ws. 

The matrix of the left Cauchy-Green tensor dj(t, r) 
becomes 

A 2 (Q 

A 2 (r) 



K{t) 

MO 



CyU, T) 



Mr) 

MO 



(2) 




I— l/A 



and 



f/ > A 2 (*), A(r) 
7i(t, i-) = r^— + 2-— 

Air) AU) 



7 2 (J, r) = 2—- 7 + T^r 



(3) 



Mr) A 2 U) 
From (1), (2). and (3) we get 
(T(t) = crn(t) — 0-22W 



■£ 



X 2 U) X(t) 

x 2 (t) MO 



/? ('M, ,_ t ), /t 

A(t) 



where 






(dU Mr)d 
U MOW 



Thus, if we know the function h, we can predict the 
stress response to an uniaxial strain history. 2 For a 
single step stress relaxation, for which A(*) = l for 
t < and A(0 = A, a constant, for t 5= 0, eq (4) becomes 



:(\ 2 -£)//(\, *) 



sphere 



and 



//(A 



■• w, = 



MX, 0* 



«x,.)=-^- 



(5) 



(5a) 



(5b) 



-For discussion of other strain histories we refer the reader to a recent paper by Bern- 
stein [3|. 



FIGURE 1. A representative three dimensional plot 0/ <r R — H(X, t). 



In their 1963 article Bernstein, Kearsley, and Zapas 
[1] plotted isochrones of a reduced stress ov?(A, t) 
which they defined by 



<tr(K = 



(T(t) 



It is evident from (5) that o7?(A, t) is identical with 
//(A, t). A typical three dimensional plot of H versus 
log time and 1/A is the surface shown in figure 1. We 
see, then, that with the stress histories of single step 
uniaxial strains we can measure //(A, t) and predict 
the stress history of any uniaxial strain history. 



3. Experimental Procedure 

All our experiments in this work deal with uniaxial 
deformations on strips with rectangular cross sections. 
These strips measured about 17 cm long, 0.88 cm wide, 
and 0.21 cm thick. The data were obtained on polyiso- 
butylene VISTANEX L-100 samples and on poly- 
vinyl chloride samples (plasticized with 50 percent by 
weight of tricresol phosphate). The stress-relaxation 
experiments were conducted on a relaxometer whose 
description has been given in a previous paper [1]. 
For the creep and recovery experiment we first allowed 
the sample to creep for a predetermined amount of 
time by attaching a known weight to the bottom clamp 
of the sample. We then cut the sample at a place 
very near the bottom clamp in order to assure an 
almost instantaneous release of force on the sample 
and measured the recovery. The constant rate of 
strain measurements were obtained on an 'Tnstron" 
machine. 
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For creep, recovery and constant rate of strain 
measurements, the motions of two marks on the 
sample, initially 10 cm apart, were obtained with the 
aid of two cathetometers. The reference marks were 
placed far enough away from the clamps to insure 
that end effects would not disturb the measurements. 
During the experiments two persons, one at each 
cathetometer, took readings and a third person re- 
corded the times. Through the use of graphical 
interpolation, we were able to ascertain accurately 
the deformation of the sample. 

4. Results and Discussion 

It has been reported previously by Bernstein, 
Kearsley, and Zapas [1] that the reduced stress for 
simple extension measurements can be represented 
very well by the following equation 



//(X, t) = (T R (k, t) = a(t)(X 2 -l)- 



--(3(t) + y(t). (6) 

A 



Although we may obtain this equation by assuming a 
particular form of the potential function £/, we shall 
regard eq (6) as an empirical equation which describes 
the surface H{\, t). Thus, for the sake of convenience, 
without committing ourselves to any particular form 
of f/, 3 we shall henceforth present and calculate our 
data with the three material properties a(t), (3{t), and 
y{t) of eq (6). Values of these functions are given in 
tables 1 and 2. In order to describe the data on plas- 



TABLE 1. Polyisobutylene L-100 at 25 °C 





Time 


(3(t) x 10 « 


y(t)x 10 s 


min 


Dynes/cm* 


Dyneslcm 3 


2 


1.76 


6.3 


4 


1.68 


5.1 


o 


1.63 


4.45 


10 


1.50 


3.75 


20 


1.29 


3.00 


40 


1.064 


2.38 


60 


0.940 


2.1 


100 


.770 


1.75 


120 


.70 


1 .65 



Table 2. Plasticized polyvinyl chloride at 27 °C 





Time 


ait) X 10~ 7 


0(f) X10" 7 


-y(f)Xl0- 7 


min 


Dynesfcm 2 


Dynes/cm 3 


Dynes/cm 2 


0.6 


2.53 


5.28 


4.36 


1.0 


2.38 


4.89 


4.01 


2.0 


2.21 


4.61 


3.64 


4.0 


2.08 


4.39 


3.54 


6.0 


2.02 


4.29 


3.45 


10.0 


1.94 


4.17 


3.34 


20.0 


1.85 


4.04 


3.23 


40.0 


1.78 


3.96 


3.13 


60.0 


1.74 


3.85 


3.07 


100.0 


1.70 


3.77 


3.01 


120.0 


1.68 


3.74 


2.99 
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ticized PVC we needed all the three material functions. 
For the polyisobutylene, however, only two were 
needed. This is in agreement with previous work [1]. 



4.1. Double Step 

Suppose that in a two step stress relaxation we take 
the first step at time t = with extension ratio X and 
the second step at t = t u with extension ratio 17, where 
the extension ratios are taken with respect to the unde- 
formed state. From eqs (4) and (5) we obtain for times 
greater than t u 



cr(t)- 



,2_i 



H{i], t) 



+ 



7)* 



XV-DL 



// 



t—t 



H Xt 



(7) 



In figures 2, 3, and 4 we show the results of two step 
stress-relaxation experiments for various values of 
t\ and the ratio of 17/X. The data are plotted taking 
as zero time the time of each step. This is done to 
emphasize the results at early times after the steps 
which are the more critical data. The black points 
represent the calculated values as obtained using 
eqs (7) and (5). The agreement is excellent. 

4.2. Triple Step 

Consider now a triple step stress relaxation with 
steps at times 0, ti, t 2 and extension ratios X, 77, £, 
respectively. For times t > t 2 we obtain from eq (4) 



Or(t)- 



H(^ty 



X 2 
(4- 



// 



> t — t] 



Hi*-**- 



H 



V 



->t 



■)}■ 



(8) 



In fii] 



5, the results of a triple step experiment 
are shown plotted against the time after each step. 
The black points represent the calculated values 
obtained by using eq (5) for the first step, eq (7) 
for the second step, and eq (8) for the third step. 
The agreement is very satisfactory. The 3-to 4-per- 
cent difference at the early times can be attributed 
to experimental difficulties. Indeed, because of these 
difficulties, we know the early time behavior for single 
step stress relaxation only to within a few percent, 
and it was these data that were used in calculating 
our predicted values. The loading time, though small 
(about 1 sec), may vary by a factor of two, and the 
manner in which the deformation is applied resembles 
a very fast ramp function rather than a step function. 
The temperature change at the early times varies with 
extension ratio. These effects could contribute signi- 
ficantly at early times after the steps. 
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FIGURE 2. Double step stress relaxation at 25 °C. 

First step with extension ratio A, = 1.602 from min to 60 min. Second step with exten- 
sion ratio A 2 = 1.880 from 60 min to 130 min. 



log t 



FIGURE 4. Double step stress relaxation at 25 °C. 

First step with extension ratio \, = 1.582 from min to 40 min. Second step with exten- 
sion ratio A.2 = 1.870 from 40 min to 165 min. 




t (min) 
FIGURE 3. Double step stress relaxation at 25 °C. 

First step with extension ratio \i = 1.272 from min to 80 min. Second step with exten- 
sion ratio k> = 1.698 from 80 min to 125 min. 



4.3. Creep and Recovery 

In a creep experiment a constant stress or a constant 
load (as in our experiment) is applied, and the resulting 
deformation is measured as a function of time. A pre- 
diction of the deformation with respect to time requires 
an inversion of eq (1) or, for the case of simple exten- 
sion, eq (4). Because this is an arduous task even 
with the aid of a computer, we elected to check our 



FIGURE 5. Triple step stress relaxation at 25 C. 

First step with extension ratio k = 1.255 from min to 20 min. Second step with exten- 
sion ratio k = 1.483 from 20 min to 60 min. Third step with extension ratio k = 1.940 from 
60 min to 125 min. 



results by putting into the right-hand side of eq (4) 
the measured deformation history up to time t and cal- 
culating the stress. 



Equation (4) can be written as 



<T(t) = 



m)- 



i 

Mt). 



H(Kt) + 
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k 2 (t) 
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(9) 
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FIGURE 6. Creep of polyisobutylene L-100 at 25 °C. 
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FIGURE 8. Creep and recovery of polyisobutylene at 25 °C. 

The arrows point to the values of the load calculated using the measured history. The The black circles indicate calculated values. White circles represent measured values, 

actual load was 86 g. 
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Figure 7. Creep of plasticized polyvinyl chloride at 27 °C. 

The arrows are used as in figure 6. The actual load was 1270 g. 

The first term of the right-hand side of this equation 
is the stress necessary to hold a material at an exten- 
sion ratio X at time t after a single step deformation 
(vide eq (5)). The second term may be taken as a 
correction for deviations of the history of the defor- 
mation from that of a single step stress-relaxation 
experiment. 

In figure 6 and figure 7 we exhibit the results of 
creep experiments on polyisobutylene and plasticized 
PVC. The line shows the value of X as a function of 
time interpolated from direct measurements. The 
arrows point to values of the load predicted by eq 
(9) and the measured history of X up to the indicated 
times. For PIB L-100 the measured load was 86 g 
and at the early times we again observe a 2-percent 
difference from the predicted load, while at the later 
times the agreement is excellent. The same is ob- 
served in figure 7 for plasticized PVC. for which the 
measured load was 1270 g. Again, this is to be ex- 
pected, since we do not know //(X. t) accurately at 
the early times. 




FIGURE 9. Constant rate of strain of polyisobutylene at 21 °C. 

The circles represent calculated points. 



During the recovery after creep, the part of the his- 
tory for which k(t) < X(r) may be considered as a rela- 
tive compression. Since the part of the surface 
//(X, t) that we measured from stress relaxation in- 
volved only extensions with X > 1, we made the as- 
sumption in our recovery calculations that eq (6) could 
be used to describe H{k. t) for values of X<1. In 
figure 8 we show the data obtained on creep and 
recovery. The black points are values of X during 
recovery, calculated from the measured values of 
X prior to recovery using eq (9) with a set equal to 
zero. 

4.4. Constant Rate of Strain 

In experiments of this type, one of the clamps hold- 
ing the sample is moved at a constant speed, while 
the other remained stationary. If there were no end 
effects due to clamping, the sample would be deformed 
at a rate proportional to the rate of separation of the 
clamps. In the experiment whose results we show in 
figure 9, we observed a constant rate of strain different 
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from what would be predicted by ignoring end effects. 
The open circles represent values of the stress cal- 
culated from eq (9). Here again, in order to predict 
the part of the curve where the motion is reversed, 
we made the same assumptions on extrapolation of 
//(A, t) to the k region as we did for the creep recovery. 
Further, these data were obtained at a lower tempera- 
ture than the stress-relaxation data. , In order to 
correct for this temperature difference, we adjusted 
the stress-relaxation data using time-temperature 
superposition appropriate to the BKZ elastic fluid 
[2]. The //(A, t) surfaces at absolute temperatures 
Tq and 7\ are related as follows 4 



The excellent agreement between the calculated 
and experimental values for several quite different 
time histories may be taken as establishing the form 
of the nonlinear dependence of stress on strain history 
given by eq (4) for elastomers in the rubbery region. 
The results reported here, being confined to simple 
extension, cannot be. taken as establishing the poten- 
tial function, U. They should be considered only as 
a test of the adequacy of the form in which time effects 
are brought into the BKZ theory. The development 
of a potential function, U, adequate for the description 
of different geometries of deformation will require 
detailed studies of different types of deformations. 



T-*p-iH(k, t) To =T-^H(K b T t) Tl 

where po and pi are the densities at To and Fi, and br 
is a time shifting factor which can either be obtained 
experimentally or estimated from semiempirical equa- 
tions. From previously reported work on polyiso- 
butylene by Leaderman [5] the shifting factor br was 
evaluated to be 1.45 for a shift from 25 °C to 21 °C. 
The experimental data and the calculated values 
shown in figure 9 are in excellent agreement. For 
creep, recovery, and constant rate of strain the cal- 
culation of each point involved graphical integration 
of the expression 



J o 



A 2 (r)" 



Mr) 
Ht). 



h \Wv'' 



dr. 



which is the second term of the right-hand side of 
eq (9). For values of r close to t we used values of 
h{k, t) extrapolated to small values of t. In most of 
our calculations the first factor in the integrand 
approached zero as r approached t, so that we expect 
the error due to extrapolation to be unimportant. 
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4 The validity of this relation for large deformations has heen tested and will be presented 
n a forthcoming publication. 
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